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Abstract We analyse spectropolarimetric observations of the penumbra of the 
NOAA AR 10953 at high spatial resolution (0.3"). The full Stokes profiles of the 
Fe I lines at 630. 1 nm and 630.2 nm have been obtained with the Solar Optical Tele- 
scope (SOT) on board the Hinode satellite. The data have been inverted by means of 
the SIR code, deriving the stratifications of temperature, line of sight velocity, and 
the components of the magnetic field vector in optical depth. In order to evaluate the 
gas pressure and to obtain an adequate geometrical height scale, the motion equation 
has been integrated for each pixel taking into account the terms of the Lorentz force. 
To establish the boundary condition, a genetic algorithm has been applied. The final 
resulting magnetic field has a divergence compatible with inside its uncertainties. 
First analyses of the correlation of the Wilson depression with velocity, tempera- 
ture, magnetic field strength, and field inclination strongly support the uncombed 
penumbral model proposed by Solanki & Montavon IT2l . 



1 Introduction 

The filamentary structure of sunspot penumbrae still harbours many open questions: 
e.g. the source of the Evershed flow, the origin of dark cored penumbral filaments 
discovered by Scharmer et al. (1 Q, as well as the structure of penumbral filaments, 
which could be elevated magnetic features (uncombed penumbral model by Solanki 
& Montavon lfl2l . Borrero et al. (f2), 0), convective penetrations (gappy penumbral 
model by Spruit & Scharmer iflOl ). or a combination of both scenarios, as suggested 
by MHD simulations by Rempel et al. Sophisticated inversion methods give 
information about the stratification of physical quantities like temperature T, mag- 
netic field strength B, field inclination y, field azimuth 0, and line of sight velocity 
Vios versus optical depth, but do not provide information about geometrical heights 
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(the Wilson depression can not be directly obtained by the inversions). Carroll & 
Kopf [I] obtain stratifications of physical quantities in geometrical height apply- 
ing a neuronal network inversion technique based on MHD simulations of quiet 
sun. However, only information included in the MHD simulation can be retrieved. 
The establishment of a geometric height scale is important for the determination of 
the electrical current vector J (crucial for the determination of ohmic energy dis- 
sipation and important for the extrapolation of the magnetic field vector from the 
photosphere to the chromosphere), for the modelling of magnetic features in 3 di- 
mensions, and for a demonstration of the reliability of MHD simulations. 



2 Observations 

The active region AR 10953 was mapped at an heliocentric angle = 10° using 
the spectropolarimeter of the Solar Optical Telescope on board the Hinode space- 
craft (Lites et al. @; Kosugi et al. 0]) on I st of May 2007, between 10:46 am and 
12:25 am UT. The region was scanned in a thousand of steps, with a step width of 
0.148" and a slit width of 0.158", recording the full Stokes vector of the pair of 
the neutral iron lines at 630 nm with a spectral sampling of 21.53 mA. The spatial 
resolution was ~ 0.32". The integration time was 4.8 s, resulting in an approximate 
noise level of 1.2 x 10~ 3 . The wavelength calibration has been done, assuming that 
the average umbral profile does not exhibit velocities. In the further study we cen- 
tre on a uniform penumbral part of a large sunspot with radial aligned filaments, 
presented in Fig.Q] 
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Fig. 1 Continuum intensity 
at 630.2 nm of a fragment of 
the active region AR 10953 
observed by SOT on board the 
Hinode spacecraft. The white 
rectangle marks the area under 
study in the present paper. 
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3 Inversion procedure 
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To derive the physical parameters of the solar atmosphere as a function of continuum 
optical depth, i.e. temperature T{x), magnetic field strength B{x), field inclination 
y(f), field azimuth <p (t), and line of sight velocity ViosIt), the SIR code (Ruiz Cobo 
& del Toro Iniesta |8|) has been applied on the spectropolarimetric data set. The 
values of those parameters are retrieved at a number of optical depth points called 
nodes. Considering the high-spatial resolution of Hinode observations and assum- 
ing that penumbral structures are horizontally resolved, an inversion with only one 
component has been performed, allowing 5 nodes in T(t), 3 nodes in B(z) and 
Vios(t), and 2 nodes in y(f) and 0(t). We have neither considered microturbulent 
velocities nor stray light contamination. The final synthesised profiles have been 
convolved with a macroturbulence velocity V mac as an additional free parameter of 
the inversion. 



4 Determination of a geometrical height scale 

The inversion method delivers for each pixel the stratification of an atmospheric 
model vs. continuum optical depth, i.e. we obtain R(x,y, t), T(x,y, t), etc. . . A geo- 
metrical height scale z(x,y, t) can be derived by the integration of 

di=-Kpdz. (1) 

For this integration three ingredients are needed: K (continuum absorption coeffi- 
cient per gram), p (mass density) and the boundary condition Zu/=z(t=1) (Wilson 
depression). 

The first ingredient, k, is evaluated by SIR from temperatures and abundances. 
It obviously depends also on gas pressures. However, it can be proven, that this 
dependence is so small that accurate K values can be obtained from very bad P g 
values. In fact, SIR obtains the gas pressure stratification by the integration of the 
hydrostatic equilibrium equation in optical depth scale. 

The pressure stratification is used to calculate the second ingredient, the density, 
using the state equation of an ideal gas and considering partial ionisation. 

If we set Zw = at all pixels, a geometric height scale z(x,y, t) can be constructed 
after integration of Eq.Q] B(x,y,z) maps are obtained by interpolation of B(x,y, t) 
maps resulting from the SIR inversion. Obviously, such R(x,y,z) maps have a di- 
vergence different from zero. Furthermore, the models are not in mechanical equi- 
librium. An optimum choice of Zw{x,y) would minimise at a given height both, the 
divergence of the magnetic field VB and the error in the motion equation. 

Neglecting viscosity, the motion equation can be written down as 



F = JxB + pg-VP ? 



(2) 



4 



K. G. Puschmann, B. Ruiz Cobo, V. Martinez Pillet 



If we neglect accelerations, F must be zero. To ensure the physical meaning of the 
solution, also VB must be zero. 
We define a merit function as 

X 2 = £ wi ( F ? + F y + F?) + w 2 (VBVB) , (3) 

pixels 

where w\, m>2 are coefficients introduced to adequately weight the contribution of 
both, VB and F. By introducing vertical displacements of the atmospheric models at 
each pixel we try to minimise the merit function (Eq.0 evaluated at a height level of 
200 km, where the uncertainties of the inverted magnetic field are minimal. We use 
a genetic algorithm changing a chromosome (array D z (x,y) containing the displace- 
ments of the atmospheric models at each pixel) until the merit function is minimised. 
The genetic algorithm has been kindly provided by Paez Mafia (software engineer 
at the Instituto de Astroffsica de Canarias). Owing to the uncertainties of the inver- 
sion method, the obtained stratification of the magnetic field is neither solenoidal 
nor satisfying the motion equation. Thus, for the evaluation of % 2 we use slightly 
modified values B'(x,y,z) =B(x,y,z) + NB(x,y), Y(x : y : z) = Y(x,y,z)+N r (x,y), and 
(j)'(x,y,z) = (j)(x,y,z) +N$(x,y), with absolute values of NB(x,y), Ny(x,y), N(j,(x,y) 
smaller than the error uncertainties of the respective parameters. The resulting syn- 
thesised Stokes profiles, considering B', / and (j)', are still compatible with the ob- 
served Stokes profiles. In summary, the solution found by the genetic algorithm is 
constituted by the optimum values of D z (x,y), NB(x,y), Ny(x,y), and N<j,(x,y). The 
best solution is reached setting wiand W2 such that both addends in Eq. [5] contribute 
in an equal manner. In each realisation the code produces slightly different results 
with a Gaussian distribution around a mean value at each pixel. Therefore, we adopt 
as final solution the average of 20 individual realisations. 

After introducing D z (x,y), Ns{x,y), N r (x,y) and (x,y) and interpolating to a 
common z- scale we can suppose, that the layer at 200 km approximately satisfies 
both, VB and the motion equation, although the pressure stratification for each pixel 
continues being the HE one. We can obtain a more accurate P g - stratification by 
the integration of the z - component of Eq. [2] However, when the P g - stratification 
changes, the z - scale is also modified. Therefore it is easier to integrate this equation 
in terms of optical depth. Taking into account that both, absorption coefficient K and 
mean molecular weight jj. hardly depend on gas pressure, we can rewrite the vertical 
component of Eq. |2]in terms of T, setting F = 0: 

/iPn dP„ iiPg ^ 
K^—^ — -=,?— - JxB - . (4) 
RT dx RT ' lz ' 

After the integration of this equation, e.g. by Runge-Kutta, we obtain p and sub- 
sequently the new z - scale from Eq. Q] This procedure has to be iterated due to the 
slight modification of B and J values. The convergence is reached after just 2 it- 
erations. The models at all pixels are interpolated to a common global z - scale. In 
Fig. |2] we present the histograms of VB at four different layers. 
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Fig. 2 Histograms of the divergence of the magnetic field vector VB at four different height layers 
before (red lines) and after (black lines) the application of the genetic algorithm. Binsize of 0.1. 



5 Results and Discussion 



Using a genetic algorithm we have evaluated the optimum boundary condition for 
the integration of the vertical component of the motion equation taking into account 
Lorentz forces. The resulting atmospheric model, interpolated to a common z - scale, 
shows very small VB values in a height range from 50 km to 200 km. This demon- 
strates the consistence of the results, although the genetic algorithm just minimises 
VB at a height of 200 km. Having at our disposal a common geometrical height 
scale we can evaluate e.g. electric currents, the Wilson depression, and in general 
the three dimensional structure of penumbral features. All these topics will be ad- 
dressed in subsequent papers. Just at a first glance we present in Fig.[3]scatter plots 
of several physical quantities at given geometrical heights vs. the Wilson depression 
(Zw). We observe a strong correlation between temperature T andZ^: in places with 
higher T, log T = is displaced to upper layers due to the strong dependence of opac- 
ity with temperature. Note beside that the z = 200km layer is nearly isotherm. The 
magnetic field strength B shows a very weak correlation with Zw- higher Zw tends 
to be related to weaker magnetic fields. Nevertheless, we find clear trends between 
field inclination y and the Wilson depression, at all layers. The same behaviour is 
observed in case of line of sight velocity Vi os . In Fig. [3] we have plotted in red colour 
points corresponding to pixels harbouring significant values of the vertical compo- 
nent of the Evershed flow (Vi os < -0.2kms _1 ). Focusing on the distribution of red 
points throughout the panels we can conclude that the Evershed flow corresponds to 
areas with increased Wilson depression, hotter temperatures, and weaker and more 
horizontal magnetic fields. All these properties support the uncombed penumbral 
model proposed by Solanki & Montavon ifPH (see also Ruiz Cobo & Bellot Rubio 
0). 
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Fig. 3 Scatter plots of temperature T, magnetic field strength B, field inclination y, and line of 
sight velocity Vi os vs. Wilson depression; at z = 0km (upper panels), z = 100km (middle panels), 
and z = 200km (bottom panels). Points with velocities at s = 0km smaller than -0.2kms _1 are 
presented in all panels by red colour. 
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